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ABSTRACT 

We report the detection of the 6.4 keV Iron Ka emission Hne in two infrared-luminous, massive, star-forming 
BzK galaxies at z = 2.578 and z = 2.90 in the CDE-S. The Chandra 4 Ms spectra of 5^^4892 and 5^^8608 
show a reflection dominated continuum with strong Iron lines, with rest-frame equivalent widths EW~2.3 keV 
and 1.2 keV, respectively, demonstrating Compton thick obscuration of the central AGN. For 5z^8608 the 
line identification closely matches the existing photometric redshift derived from the stellar emission. We use 
the observed luminosities of the Iron Ka line, of the rest-frame mid-IR continuum and of the UV rest-frame 
narrow emission lines to infer intrinsic L2-10 keV ^ 10^"^ erg s~^ , about 1 .0-2.5 dex larger than the observed ones, 
hence confirming the presence of an absorber with A^h > 10^^ cm~^. The two BzK galaxies have stellar masses 
of 5 X 10^^ M0 and, based on VLA 1.4 GHz and submm 870/im observations, they appear to host vigorous 
starburst activity with SFR ~ 300-700 M© yr~^ that is also optically thick. We estimate that the AGN might also 
conceivably account for an important fraction of the bolometric far-IR emission of the galaxies. The implied 
volume density of Compton thick (CT) AGN with L2-10 keV > 10"^"^ erg s~^ is in agreement with predictions 
from X-ray background synthesis models. These sources provide one of the first clearcut observations of the 
long-sought phase of simultaneous, heavily obscured quasar and star formation activity, predicted by models 
of massive galaxy evolution at high redshifts. 

Subject headings: galaxies: general — galaxies : active 



1. INTRODUCTION 

Highly obscured, deeply embedded AGN might contribute 
significantly to the total accretion power in the Universe (Mar- 
coni et al. 2004) and are required to account for the spectrum 
of the X-ray background (XRB, GilH et al. 2007). In partic- 
ular, AGN obscured in the X-ray band by column densities 
larger than N// = 10^"^ cm~^ (Compton thick AGN) represent 
20-25% of the AGN detected by Integral and SWIFT/BAT in 
the local universe (Malizia et al. 2009; Burlon et al. 2010). 

The identification of highly obscured AGN is particularly 
challenging since, because of the high column densities, the 
vast majority of their X-ray emission below rest frame 10 keV 
is absorbed, so that they are largely missed even in the deep- 
est hard X-ray surveys available today. At high redshift the 
quest for such objects is mainly pursued by indirect evidence, 
by selecting galaxies with high ratios of mid-infrared to op- 
tical, or mid-infrared to X-ray fluxes, and by deriving aver- 
age X-ray properties by stacking techniques (e.g., Daddi et al. 
2007b; D07 hereafter; Fiore et al. 2008, 2009; Alexander et 
al. 2005; 2008; Treister et al. 2010, Georgantopoulos et al. 
2009, Eckart et al. 2010, Donley et al. 2010). These studies 
have suggested that a substantial fraction of massive galaxies 
at z > 1 host highly obscured, intrinsically luminous AGN. 

Only a few tens of secure Compton thick AGN are known 
so far in the local Universe (Comastri et a. 2004, Dalla Ceca 
et al. 2008 and references therein), and only for a handful of 
them a quasar-like intrinsic hard X-ray luminosity has been 



|chiara .feruglio @ cea.fr I 

^aboratoire AIM, CEA/DSM - CNRS - Universite Paris Diderot, 
Irfu/Service d'Astrophysique, CEA Saclay, Orme des Merisiers, 91191 
Gif-sur-Yvette Cedex, France 

^ INAF - Osservatorio Astronomico di Roma, via Frascati 33, 00040 
Monteporzio Catone, Italy 

^ Department of Physics, Durham University, South Road, DHl 3LE, 
UK 

^ Universita' La Sapienza, Piazzale Aldo Moro 5, 00185 Roma, Italy 



inferred (L2-10 keV > lO'^'^erg s~^; Braito et al. 2004, Picon- 
celli et al. 2010). Their X-ray spectra are characterized by 
the presence of a strong Iron Ka 6.4 keV fluorescent emis- 
sion line, with large equivalent width EW > 1 keV, on a flat 
reflection-dominated continuum (Matt et al. 2000). Not much 
evidence currently exists in the distant Universe of AGNs 
showing prominent, high EW Iron Ka line: an IRAS selected 
hyperluminous galaxy at z = 0.93 (Iwasawa et al. 2005) and 
two possible detections in the CDFS reported by Norman et 
al. (2002; a type 2 QSO at z = 3.7) and by Tozzi et al. (2006; a 
galaxy at z = 1 .53); both confirmed by Comastri et al. (201 1). 

In the following we present the significant Chandra de- 
tection of the Iron Ka emission line in two massive BzK 
galaxies at redshift z= 2.578 and z= 2.90 ± 0.10 in GOODS- 
South, using the recently acquired 4 Ms dataset. We used 
for this purpose the entire set of 52 observations for a to- 
tal exposure time of 4 Ms, available in the Chandra Data 
Archive (http://cxc.harvard.edu/cda/Contrib/CDFS.html). We 
adopt a ACDM cosmology (Hq = 70 km s"^ Mpc"^ Qm=03; 
= 0.7) and a Chabrier stellar initial mass function (IMF). 

2. X-RAY PHOTOMETRIC AND SPECTRAL ANALYSIS 

We have been investigating the properties of BzK -selected 
(Daddi et al. 2004), massive galaxies at redshift 1 < z < 3 
with faint or no X-ray emission, with the aim of constraining 
their obscured AGN activity. We have focused in particular 
on objects showing also a power-law emission in the near- 
IR rest frame, an independent indication for the presence of a 
luminous, obscured AGN. A large statistical sample of 750 X- 
ray undetected galaxies has been assembled in the COSMOS 
field, and their X-ray properties investigated through stacking. 
As a control sample, we also studied in detail a smaller sam- 
ple of 26 similarly selected galaxies in GOODS-South, with 
less than 200 net counts in the 0.5-7 keV band in the CDFS 
4 Ms data that would remain thus undetected in the shallower 
COSMOS Chandra observations. While the general results 
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Fig. 1 . — The spectra of BzK4^92 (upper panels) and BzK^60^ (lower panels), fitted with a pure reflection model plus an Iron Ka 6.4 keV line. The spectra 
are shown both in observed counts units (left panels) and in physical units (keV^; right panels). Each bin corresponds to a 2a measurement. The red dotted line 
represents the average background extracted from the whole field, excluding sources (Fiore et al. 2011, in preparation). 



of these studies will be reported elsewhere (Feruglio et al. 
2011, in preparation), we concentrate here on two remark- 
able sources in our GOODS-S control sample, for which the 
available Chandra data allow the unambiguous detection of 
the Iron Ka emission. 

Both of the BzK galaxies under exam were selected from 
the D07 sample. BzK4^92 is a source with spectroscopic red- 
shift z=2.578 (Szokoly et al. 2004, Vanzella et al. 2006; based 
on multiple narrow emission lines), detected in the X-rays al- 
ready in the 1 Ms Chandra data (Giacconi et al. 2001) as 
well as in the 2 Ms data (Luo et al. 2010). In the 4 Ms 
dataset it has 75 net counts in the 0.5-7 keV band and hard- 
ness ratio of 0.58 (the hardness ratio is defined as HR=hard- 
soft/hard+soft, where S and H are the soft and hard band net 
counts detected by Chandra). BzK^60^ has a photometric 
redshift Zphot = 2.94 (from D07; consistent with Zphot = 2.88 
from Santini et al. 2009). VLT+VIMOS spectroscopic obser- 
vations (Popesso et al. 2009) show that there is no continuum 
emission in the blue, consistent with 3, but failed to yield 
a reliable redshift estimate. We tentatively detected in the VI- 
MOS spectrum a low signal to noise feature at 6027 A (a re- 
gion clear from OH lines) that, if identified with CIV, would 
give z= 2.88. We notice that the photometric redshift for 
5z^8608 is reliable as the galaxy fulfills the U-band dropout 
criteria, corresponding to a Lyman continuum break (being 
undetected in the deep U-band imaging of Nonino et al. 2010; 
Fig. [3]). This source is not listed in the Chandra 2 Ms catalog 



(Luo et al. 2010), likely because of the presence of another 
nearby X-ray source (#193). In the 4 Ms data BzK8608 has 
92 counts in the total band. It is not detected the 0.3-1.5 band, 
and its hardness ratio is HR>0.54 (3a). 

We extracted the Chandra spectra using a circular region 
of 1.5'' and 3" radius for BzK4892 and BzK8608 (off-axis 
~7.40, respectively, estimating the background in an annulus 
around the source position. We experimented with different 
extraction radii, centering, and background estimates, taking 
care in all cases of excluding nearby detected sources (partic- 
ularly important for BzK8608), finding that the features found 
in the spectra were robust to the detail extraction method. The 
spectra for both galaxies show a flat continuum and promi- 
nent emission lines peaking at ~ 1.78 keV for BzK4^92 and 
- 1.64 keV for 5z^8608 (Fig. 1). 

We consider two spectral models to fit the data in XSPEC: 
i) a reflection-dominated model, pexrav + zgauss, and ii) 
a transmission (power-law) model, wahs^(pow+zgauss). A 
narrow gaussian component is used to model the emission 
lines. Given the limited photon statistic, the slope of the 
primary continuum has been fixed to F = 1.8 in both scenar- 
ios. Based on C-statistics, the emission line is detected with 
a significance of >4cr (AC = 17.5) for BzK4'^92 and > 3cr 
(AC = 7.5) for 5z^8608 (see also Fig. [3 for the reflection- 
dominated model. In the case of a transmission model the sig- 
nificance of the line remains above 4cf for BzK4^92 and just 
below 3a for 5z^8608. To further evaluate the significance 
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Fig. 2. — Upper panels: Chandra images of BzK4^92 and BzK^60^ in three energy bands. The middle images are centered at the energy of the Iron K line. 
The circles show the extraction regions. Lower panels: confidence levels (1,2 and 3 a) obtained from the fit of the Iron line. The vertical blue lines show the 
optical spectroscopic/photometric redshift constraints. 



of the observed lines we also performed 10^ simulations of 
the source continuum and background for both sources. The 
probability that the lines might be due to chance fluctuations 
are 3x10"^ and 2x10"^, respectively. The spectral fits are 
summarized in Table 1 . 

Based on the spectroscopic and photometric redshifts, we 
identify the lines as Iron Ka for both sources. In particular, 
we derive a spectroscopic redshift for 5^^8608 from the X- 
ray spectrum, obtaining z = 2.88 ± 0.10. Figure 2 shows the 
1-2-3 a confidence level contours, showing the redshift fit and 
the significance of the Iron line detection. 

The rest frame equivalent widths of the Iron lines are 
EW= 23tli keV for BzK4^92 and EW= 1.2 ±0.4 keV for 
5z^8608 (for a reflection model, see Table 1), implying ab- 
sorbing columns exceeding 10^^ cm~^. 5^^8608 shows an 
additional low-significance excess at 1.8 keV that would cor- 
respond to the energy of Iron K/3 at z = 2.9, if real. The de- 
tection of Iron Ka at 6.4 keV excludes that the high hardness 
ratios in these sources might be due to high mass X-ray bi- 
naries, as suggested by Donley et al. (2008), due to spectral 
incompatibility (Persic & Raphaeli 2002). 

The fit with the transmission model gives N// = 7 ± 2 x 10^^ 
cm-2 for BzK4^92 and Tl x 10^^ cm'^ for 5z^8608, and 
is statistically indistinguishable from the fit with a reflection 
model. However, these Nh values are inconsistent with the 
measured EW>1 keV. Indeed, for N// ~ 6 x 10^^ cm~^, a 
EW~300 eV is predicted (Ikeda et al. 2009) and confirmed 
by observations (Guainazzi et al. 2005, Fukazawa et al. 2010). 
EW> 1 keV are the hallmark of column densities larger than 
10^"^ cm~. Therefore, in the following, we adopt as best fit the 
results from the reflection dominated model for both sources 
(Tab. El). 

Figure 2 shows the Chandra images of BzK4^92 and 
BzK^60^ in three energy bands, including one centered on 
the Iron lines. We note that 5^^8608 is not detected in the 
soft band. 

3. DISCUSSION 



3.1. Intrinsic luminosities and obscuration 
The observed (obscured) hard X-ray luminosities are 
^2-10 kev = 10"^^-^ and lO^^-^ erg s'^ for BzK4^92 and 
5z^8608, respectively. Given the large A^h implied by the 
detection of high EW Iron Ka lines, the X-ray emission is 
optically thick and we cannot directly estimate the intrinsic 
^2-10 kev using the X-ray data (e.g.. La Massa et al. 2011). 
In order to estimate the intrinsic luminosities of the AGN we 
used a variety of independent methods. 

Following Iwasawa et al. (2005), a lower limit to the X- 
ray luminosity can be obtained in the expectation that in 
typical conditions the luminosity of the Iron Ka line is at 
most 3% of the 2-10 keV luminosity. For BzK4S92 = 
1.25 X 10"^^ erg s~^ would imply a intrinsic 2-10 keV lumi- 
nosity of 4x 10"^^ erg s~^. The relation L/^q,/L2-io keV of Lev- 
enson et al. (2006) implies an higher luminosity, L2-10 keV ^ 
j^Q44.8 gj,g g-i alternative estimate is obtained from the 
6.0/im luminosity. The 24 /am flux of 591/iJy implies 
U.o/im = 10"^^-^ erg s-^ (the mid-IR SED of Mrk 231 is used 
for the small K-correction). While 24/im at z = 2.578 might 
be contaminated by PAH features, we find similar results if 
we use the 16/im flux of 209/iJy (Teplitz et al. 2010) for 
which no contamination is expected (4. 5 /im rest- frame). This 
would suggest L2-10 keV = 10"^^^^-^ erg s~^ using the Lutz et al. 
(2004) correlation and a lower L2-10 keV = lo"^"^ "^^^-^ erg s~^ 
using Lanzuisi et al. (2009). A third estimate can be derived 
from the luminosity of the UV-optical emission lines. The 
VLT/FORSl-2 optical spectrum of BzK4S92 shows promi- 
nent Lya, CIV, Hell and GUI] narrow emission lines (Szokoly 
et al. 2004, Vanzella et al. 2006). Using the observed ratios 
between the line and the intrinsic 2-10 keV emission from 
Mulchaey et al. (1994), we infer L2-10 keV = lO^^-^iO-S erg s"^ 
Using the Netzer et al. (2006) conversion for UV-optical line 
luminosities would give larger X-ray luminosities by ~0.5 
dex, or L2-iokeV = 10^^^^-^ erg s~^ All in all, it appears 
that the 2-10 keV luminosity of BzK4^92 is in the QSO 
range, well in excess of 10^^ erg s~\ and more likely of 
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TABLE 1 
Spectral fit parameters 



Source ID 


Model 


Fe Ka flux 
10"^^ erg/cm^s 




EW 
keV 


Cstat/d.o.f. 


BzK4892 
BzK4892 


Abs 
Refl 


1.6zb0.5 
2.0dz0.6 




2.1=b0.5 

^•-^-0.6 


307/441 

297/442 


BzK8608 
BzK8608 


Abs 
Refl 


3.6zb0.5 
3.4=b0.8 


7.0!|g 


l.OzbO.3 
1.2zb0.4 


411/440 
418/442 



^q44.5 45 gj,g g 1 jj^g implied column density is at the level of 

A^H - 1024-5-25 cm-l 

For 5zi^8608, we measure L(Ka)= 2.5 x lO^^ erg s'^ . This 
converts to L2-10 keV ~ lO"^"^"^ erg s~^ (Iwasawa et al. 2005), or 
~ 10^5 erg s~^ adopting Levenson et al. (2006). This galaxy 
is about 10 times fainter in the mid-IR than 5^^4892, which 
implies a lower L6.0 fim- Using the Lutz et al. (2004) rela- 
tion, this in turn suggests L2-10 keV = 10"^"^ erg s~^. We cannot 
obtain an independent estimate from the UV-optical emission 
lines given that no lines were detected in the VIMOS spec- 
trum. The intrinsic 2-10 keV luminosity of BzK^60^ is likely 
^ 10^4 erg s~^. The implied column density is likely in the 
range A^h 1-5 x lO^^ cm"^, hence this object might be a bor- 
derline CT AGN. We notice that both sources have estimated 
intrinsic 2-10 keV luminosities in the QSO regime. 



3.2. Bolometric luminosities and star formation rates 

The ratio of the bolometric to the 2-10 keV luminosity, 
Lz,o//L2-io kev» is typically of the order 30-50 for quasars 
(Marconi et al. 2004). Assuming such a bolometric correc- 
tion, we estimate U^/ - lO^^iO-S ^-i BzK4S92 and - 
j^Q45.5±o.5 gj,g g-1 5^^8608. These very high luminosities, 
if expressed in solar units, correspond to Lhoi ^ 
and ~ 10^^^^^ ^L©. When powered by star formation, simi- 

TABLE2 
Galaxy properties 



BzK4S92 BzKSGOS 



Ra 


3:32:35.71 


3:32:20.95 


Dec 


-27:49:16.1 


-27:55:46.3 


Redshift 


2.578 


2.90zb0.10 


HR 


0.58 


>0.54 (3a) 


SFRuv [M© yr"^] 


70 


70 


SFRRadio[M0 yr"^] 


1100 


300 


L(Ka) [ergs-^] 




1042.4 


L2-10 kev, abs [erg s"^] 




1042.8 


L2-iokev, K« [ergs~^] 


-j^Q43.6-44.8 


1Q43.9-45 


L2-10 kev, 6/2m [erg s"^] 


-j^q44.4-45 


1^44 


L2-IO keV, UV lines [crg S~^] 


^q44-45.5 




i [AB mag] 


24.75 


24.72 


IRAC 3.6 /im [AB mag] 


20.81 


21.29 


IRAC 4.5 /im [AB mag] 


20.54 


21.16 


IRAC 5.8 /im [AB mag] 


20.1 


20.84 


IRAC 8.0 /im [AB mag] 


19.64 


20.65 


MIPS 24 //m [mJy] 


0.591zb0.07 


0.045zb0.003 


MIPS 70 /im [mJy] 


3.3zbl.8 




850 /im [mJy] 


3.3zbl.l 




1.4 GHz [/iJy] 


80zbl4 


17.5zb6.5 



lar luminosities require star formation rates (SFR) in the host 
galaxies at the level of 100-1000 Mq yr~^. It is interesting 
thus to compare these estimates to the inferred SFR for the 
two BzK galaxies, in order to compare the relative AGN and 
SFR contributions. 

From the stellar SED (UV to the near-IR rest frame; Fig. 3), 
we estimate SFR ~ 70 Mq yr~^ for 5^^4892, corrected for 
dust reddening based on the observed UV slope (see D07 for 
more details). It is quite possible though that also the stellar 
UV emission might be optically thick. If we interpret its mid- 
IR emission (591 /iJy at 24/im) as due to star formation this 
would imply a whopping SFR~ 10"^ Mq yr~\ demonstrating 
that the mid-IR is likely completely dominated by the AGN 
(this object is among the most extreme mid-IR excess galax- 
ies in the sample of D07). The SED fit with a star forming 
galaxy template (Fig. 3) shows an excess emission already 
at 5.8 /im (observed frame) probably due to the AGN con- 
tribution. The SED in the IRAC bands is indeed showing a 
steep power law, extending all the way to 24 /im. Due to its 
faint optical counterpart (i~25 AB), BzK4^92 has high mid- 
infrared to optical flux ratio, F(24)/F(R)~2000 and therefore 
is also classified as a dust-obscured galaxy (DOG, Dey et al. 
2008). It is an extreme object also in the Fiore et al. (2008) 
sample. Alonso-Herrero et al. (2006) reports a 70 /im flux 
of 3.3 ± 1.8 mJy for this galaxy, which is also likely affected 
by the AGN emission. On the other hand, the galaxy is seen 
at 1.4 GHz with a flux of 80 /iJy in the VLA data of IVliller et 
al. (2008). If due to star formation, using the radio-IR corre- 
lation, this would imply SFR~ 1 100 Mq yr~^. Inspecting the 
pubhcly available Apex-^LABOCA 870/im map of GOODS- 
S (Weiss et al. 2009), we find a 3a signal at the position of 
BzK4S92 of 3.3 mJy, which would also convert into a similar 
SFR, -500 1M0 yr-i. 

Even for the most luminous quasars it is generally found 
that the far-IR emission in the sub-mm bands is due to star 
formation. If this is the case also for BzK4^92, our results 
suggest that the galaxy is witnessing also very powerful star 
formation activity at the level of 500-1000 Mq yr"^ Only 5- 
10% of this is seen directly in the UV, implying that also the 
UV-emission from stars is optically thick, similarly to local 
ULIRGs (e.g., Goldader et al. 2002; da Cunha et al. 2010) 
and to what is expected to be found in major mergers. How- 
ever, we notice that the radio and 870/im derived SERs would 
formally imply L^oi ^ 10^^"^^^ erg s~\ comparable to what 
inferred for the obscured AGN. Hence, the AGN is likely con- 
tributing an important fraction of the total L^^/- 

The morphology of the galaxy from HST-fACS imaging is 
suggestive of a merger, showing two distinct clumps in the 
UV in the rest frame. When smoothing, after excluding the 
two bright knots, we detect significant faint low-level emis- 
sion in the summed i-\-z band, extended over a diameter of 
about 0.75'' (6 kpc). Overall, this is consistent with the size of 
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Fig. 3. — Optical to mid-infrared SEDs of BzK4^92 (left panel, black symbols) and BzK^60^ (right panel). The blue continuous line shows the fitted model 
SED of a pure star forming galaxy. The photometry from the GOODS datasets is taken from D07. HST+ACS color images are also shown as thumbnails. 



a massive galaxy at 2, and we cannot exclude that the two 
UV knots might be just luminous HII regions inside a big disk 
galaxy. 

The photometric information is of lower quality for 
BzK^60^ due to its overall faintness. The SFR inferred from 
the UV is also 70 M^/yr. This object is also a mid-IR ex- 
cess source in D07 sample (by a less extreme factor of 6), 
with S24 = 45 /iJy ( -1/10 of 5^^4892), although it is not a 
DOG. The SED (Fig. 3) might also be consistent with pure 
stellar emission, but a deviation from the star formation tem- 
plate is observed at the 8 /im band. This galaxy is not detected 
at 70 /im and 870 /im, suggesting a lower AGN luminosity 
and/or SFR. Inspecting the VLA 1.4 GHz data (Miller et al. 
2008) we detect a faint 3a source at its position with a flux 
of 17.5 ± 6.5 /iJy. This corresponds to a luminosity 3 x 10^^ 
L0 and SFR— 300 yr~^. If the radio signal is real and 
not due to AGN, also in this case the stellar emission appears 
to be heavily obscured, i.e. optically thick in the UV, as ex- 
pected for very active starbursts and mergers. The HST+ACS 
imaging of this galaxy (Fig. 3) is not particularly telling. 

For both galaxies we derive similar estimates of the stellar 
masses, at the level of 5 x IO^^Mq. Our results thus support 
the picture that the luminous, Compton thick AGN that we 
discovered are hosted by fairly massive galaxies at z ~ 2.5-3, 
that at the same time also host vigorous star formation activ- 
ity heavily obscured by dust. This is relevant in the AGN-host 
galaxy co-evolution scenario, in which a phase of rapid, heav- 
ily obscured BH growth accompanied by intense obscured 
star-formation is predicted (Silk & Rees 1998, Fabian 1999, 
Granato et al. 2004). 

3.3. Implications for Compton thick nuclear activity at high 

redshift 

We derive a crude value of the volume density of Compton 
thick AGN with high EW Iron Ka emission, based on our 2 
detections. We conservatively use the full redshift range ex- 
plored by the BzK selection (z=l .4-3), finding a space density 
of 3 X 10~^ Mpc~^. This is consistent with the predictions of 
the Gilli et al. (2007) X-ray background synthesis model, for 



Compton thick AGN with L2-10 kev > 10^^ erg s~^. However, 
we notice that our sampling of L2-10 keV > 10^^ erg s~^ AGN 
with Compton thick absorption at 1 .4 < z < 3 might be sub- 
stantially incomplete due to selection effects, and the lumi- 
nosity of at least one of the BzK galaxies in this paper might 
be substantially higher than the adopted 10"^"^ erg s~^ limit. 

We emphasize that the sources presented here have X-ray 
spectra that are the high-luminosity analogous of those of lo- 
cal prototype Compton thick AGN. This discovery represents 
one of the first clear-cut evidence for this class of objects at 
high redshift. Our results confirm that indeed there is a pop- 
ulation of Compton thick QSO among massive, star-forming, 
dust obscured galaxies, as suggested by several studies (D07, 
Fiore et al. 2008, 2009, Treister et al. 2010, Alexander et al. 
2008, Lanzuisi et al. 2009), and as predicted by models of 
merger-driven AGN/galaxy co-evolution (Silk & Rees 1998, 
Fabian 1999, Di Matteo et al. 2005). Indeed, these lumi- 
nous Compton thick quasars appear to be hosted by galaxies 
with optically thick UV emission, which are a relative small 
minority among BzKs (Daddi et al 2007a), and imply the 
presence of a dense starburst usually connected with merg- 
ers. This seems to agree with evolutionary scenarios, where 
mergers induce both rapid accretion onto SMBH and intense 
star-formation activity. 

The two sources presented here have high intrinsic X-ray 
and bolometric luminosity, and bright Iron Ka line, allowing 
their identification. They might represent just the tip of the 
iceberg of a much larger population of highly obscured AGN 
with somewhat lower intrinsic luminosity, whose Iron lines 
are still remaining hardly detectable even in the deepest hard 
X-ray surveys. For example, the Gilli et al. model predicts 
a ^ lOx higher space density of Compton thick sources with 
L2-10 keV > 10"^^ erg s~^. Detecting Iron lines in galaxies with 
such luminosities might be challenging with current facilities, 
even if they have similar EW as in BzK4^92 and 5z^8608. 



We thank the anonymous referee for useful comments. 
This research was supported by the ERC-StG grant UPGAL 
240039 and by the French ANR under contract ANR-08- 



6 



Femglio et al. 



JCJC-0008. EP acknowledges financial support from ASI (grant 1/088/06/0). DMA thanks the Science and Technology 

Research Council for support. 



Alexander,D.M., Bauer, F.E., Chapman, S.C. et al. 2005, ApJ, 632, 736 
Alexander,D.M., Chary, R.R., Pope, A., 2008 

Alonso-Herrero, A., Perez-Gonzalez, P. G., Alexander, D. M. et al. 2006, 
ApJ, 640, 167 

Braito V., et al., 2004, A&A, 420, 79 

Burlon, D. et al. 2010, ApJin press, arXiv : 1012.0302] 

Comastri A., Brusa M., Civano R, 2004, MNRAS, 351, L9 

Comastri, A. et al. 201 1, A&A, in press 

Da Cunha, E. et al. A&A, in press, arXiv: 1008 .20001 

Daddi, E. et al. 2007a, ApJ, 670, 156 

Daddi, E. et al. 2007b, ApJ, 670, 173 

Delia Ceca R., et al., 2008, A&A, 487, 119 

Dey, A. et al. 2008, ApJ, 667, 943 

Di Matteo, T., Springel, V., Hernquist, L. 2005, Nature, 433, 604 
Donley, J.L., Rieke, G.H., Perez-Gonzalez, RG., Barro, G., 2008, ApJ, 687, 
111 

Donley, J. L., Rieke, G.H., Alexander, D. M. et al. 2010, ApJ, 719, 1393 

Fabian, A.C. 1999, MNRAS, 308, L39 

Fiore, R, Grazian, A., Santini, R et al. 2008, ApJ, 672, 94 

Fiore, F. et al. 201 1, in preparation 

Fukazawa et al. 2010, in press (arxiv 1011.1750) 

Georgantopoulos, I. et al. 2008, A&A, 484, 671 

Georgantopoulos, I. et al. 2009, A&A, 507, 747 

Giacconi, R. et al. 2001, ApJ, 551, 624 

Gilli R., Comastri A., Hasinger G. 2007, A&A, 463, 79 

Goldader, J. D. et al. 2002, ApJ, 568, 651 

Granato, G. L., De Zotti, G. , Silva, L., et al. 2004, ApJ, 600, 580 



Ikeda, S. et al. 2009, ApJ, 691, 608 

Iwasawa, K., Crawford, C.S., Fabian, A. C. et al. 2005, MNRAS, 362, 20 
Guainazzi, M., Matt, G. and Perola, G. C. 2005, A&A, 444, 119 
Lanzuisi, G., Piconcelli, E., Fiore, F. et al. 2009, A&A, 498, 67 
Levenson, N. A., Heckman, T. M., Krolik, J. H. et al. 2006, ApJ, 648, 111 
Luo, B. et al. 2010, ApJS, 187, 560 
Lutz D. et al. 2004, A&A. 418, 465 
Malizia, A. et al. 2009, MNRAS, 399, 944 
Marconi, A. et al. 2004, MNRAS, 351, 169 
Matt, G. et al. 2000, MNRAS, 318, 173 

Miller, N. A., Fomalont, E. B., Kellerman, K. I. et al. 2008, ApJS179,114 

Mulchaey, J. S. et al. 1994, ApJ436, 586 

Netzer, H., Mainieri, V., Rosati, R et al 2006, ApJ, 652, 117 

Nonino, M. et al. ApJS, 183, 244 

Piconcelli, E. et al. 2010, ApJ, 722, 147 

Persic, M. & Raphaeli, Y. 2002, A&A, 382, 843 

Popesso, P., Dickinson, M., Nonino, M. et al. 2009, A&A, 494, 443 

Santini, R et al. 2009, A&A, 504, 751 

Silk, J. & Rees, M. J. 1998, A&A, 331, LI 

Szokoly, G. R, Bergeron, J., Hasinger, G. et al. 2004, ApJ S, 155, 271S 
Teplitz, H., Chary, R., Elbaz, D., et al., 2010, AJ in press jarXiv: 1010. 17971 
Tozzi, R , Gilli, R., Mainieri, V. et al. 2006, A&A, 451, 457 
Treister, E. et al. 2010, ApJ, 722, 238 

Vanzella, E., Cristiani, S., Dickinson, M. et al. , 2006, A&A, 454, 423 
Weiss, A. et al. 2009, ApJ, 707, 1201 



